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We study the anisotropic efTect of the Coulomb interaction on a 1/3-filling fractional quantum 
Hall system by using exact diagonalization method on small systems in torus geometry. For weak 
anisotropy the system remains to be an incompressible quantum liquid, although anisotropy mani- 
fests itself in density correlation functions and excitation spectra. When the strength of anisotropy 
increases, we find the system develops a Hall-smectic-like phase with one-dimensional charge density 
wave order and is unstable towards the one-dimensional crystal in the strong anisotropy limit. In 
all three phases of the Laughlin liquid, Hall-smectic-like, and crystal phases the ground state of the 
anisotropic Coulomb system can be well described by a family of model wavefunctions generated by 
an anisotropic projection Hamiltonian. We discuss the relevance of the results to the geometrical 
description of fractional quantum Hall states proposed by Haldane [Phys. Rev. Lett. 107, 116801 
(2011)]. 
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I. INTRODUCTION 

The fractional quantum Hall (FQH) effect at an odd 
denominator filling of u has been understood as a prop- 
erty of incompressible quantum liquid in an interacting 
two-dimensional (2D) electron system. Laughlin's trial 
wavefunctioni is the first successful theory to describe 
this many-body effect, where the interacting system is 
implicitly assumed to be isotropic. Since then most the- 
oretical works on the FQH system have followed this sim- 
ple assumption and the FQH states are considered to be 
isotropic with rotational symmetry. 

However, the real FQH systems may be anisotropic. 
One natural source for this is the anisotropic dielectric 
tensor, which in turn leads to anisotropic Coulomb inter- 
action. Other mechanisms for various anisotropic FQH 
systems have also been discussed theoretically^"^ and 
experimentally i^ For example, an anisotropic FQH state 
in a z^ = 7/3 system has been observed in experimenti^ 
In these anisotropic FQH systems, rotational symmetry 
of the consequent ground state is expected to be broken. 
In the extreme anisotropic interaction limit, where the 
Coulomb interaction can be effectively treated to be one- 
dimensional (ID), the ground state of the system will 
be a quasi-lD crystal)^ Thus, the properties of the FQH 
state in an anisotropic interaction may not always be as- 
sociated with the isotropic incompressible liquid. A com- 
prehensive investigation on the effect of the interaction 
anisotropy is called for. 

Motivated by the anisotropic transport properties ex- 
perimentally reported at the partially filled higher Lan- 
dau level (LL), trial wavefunctionsiSrJ^ have been pro- 
posed to describe the anisotropic FQH states. These 
variational wavefunctions modify the original isotropic 
Laughlin wavefunction by splitting the multiple-order ze- 
ros in the wavefunction. Very recently, Haldane- has 



constructed a family of the Laughlin states, which are 
the exact ground states of the corresponding projected 
Hamiltonians and can be parameterized according to the 
interaction anisotropy. This variational state can be 
compared numerically with the anisotropic FQH state. 
An effort to map the underlying wavefunction of this 
variational state has been reported.— In this paper we 
study FQH states in 2D electron systems with anisotropic 
Coulomb interaction and discuss the relevance of our re- 
sults with the geometric description of the FQH states.— 
The paper is organized as follows. In Section |lT] we 
introduce our model with anisotropic Coulomb interac- 
tion and set up the Hamiltonian on a torus geometry. In 
Section IIIII we discuss the properties oi v = 1/3 FQH 
states at different regimes of the interaction anisotropy 
using energy spectra, charge density and correlation func- 
tions. We also compare the anisotropic FQH state with 
variational Laughlin states using wavefunction overlap. 
Section [IVI summarizes the paper. 



II. MODEL AND NUMERICAL SETUP 

We study a 2D electron system under a perpendicular 
magnetic field B = Bz. The electron-electron Coulomb 
interaction with an in-plane biaxial dielectric tensor has 
the form 



Vcir) = 



Aney/AcX^ +y^/Ac 



(1) 



where Ac is the interaction anisotropy parameter and di- 
rections of X and y are along the two principal axes of the 
dielectric tensor. The effective mass tensor is considered 
isotropic so that non-interacting electrons move in the 
circular cyclotron orbitals. However, equipotential lines 



of the Coulomb interaction are generally elliptical with 
Ac y^ 1. In the following discussion, we choose A^. > 1 
such that X is the hard axis. For Ac < I, one simply 
swaps the easy and hard axes. At Ac ^ 1, the Coulomb 
interaction is effectively a ID repulsion along the hard 
axis. 

In our numerical calculations, we use Landau gauge 
(0, Bx) for the magnetic vector potential. Periodic 
boundary conditions for the magnetic translational op- 
erators are imposed with a quantized flux N^ through 
the rectangular unit cell Lx x Ly . The magnetic length £ 
is taken as the unit length and the energy is in units 
of e^/47re^. To reduce the size of the Hilbert space, 
we carry out our calculation at every pseudomomentum 
K — (KxTKy),^"^ where Kx (Ky) is in units of 2tt/Lx 
(27r/Ly). The magnetic field is assumed to be strong 
enough so that the spin degeneracy of the Landau levels 
is lifted.— li^ One can thus project the system Hamil- 
tonian into the valence Landau level.-^ For the lowest 
Landau level, the projected Hamiltonian has the form 



Hc = ^J2n^)^'''^'T.e'''-^'"'''''\ 



(2) 



i<j 



where the momentum q = {qx,%) takes discrete val- 
ues suitable for the lattice of the unit cell and Ri 
is the guiding center coordinate of the z-th electron. 

V^(q) ~ 1/ Jqx/Ac + Acqy is the Fourier transform of 

the Coulomb interaction. From the geometrical point of 



view, we generalize q 



lx 



ll to ql 



9 qaqt- 



I /Ac 
Ac 



where 



(3) 



is the inverse metric for the Coulomb interaction. 



III. NUMERICAL RESULTS AND DISCUSSION 

For a ly — 1/3 FQH system, we first study the low- 
lying energy spectra using exact diagonalization method. 
Here and in following subsections, the default size of the 
i^ = 1/3 system is Nc = 10 and the default shape of the 
unit cell is square unless otherwise specified. We find 
qualitatively similar results in systems with other sizes 
and/or different shapes of the unit cell. 

Figure [T] plots the excitation energy gap as a function 
of the Coulomb interaction anisotropy. In the range up to 
Ac = 100 the curve is nonmonotonic and develops sev- 
eral distinct regimes. For small interaction anisotropy 
up to Ac = 2.0, the energy gap remains nearly con- 
stant, indicating the incompressible liquid phase in the 
isotropic case (i.e., Ac = 1) is robust against weak in- 
teraction anisotropy. When the interaction anisotropy 
further increases, the energy gap decreases to a min- 
imum at around Ac = 8.0. For larger Ac the energy 
gap increases with the interaction anisotropy to a maxi- 
mum at around Ac = 33.0. The finite-size scaling shown 
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FIG. 1: (Color online) Excitation energy gap versus interac- 
tion anisotropy for A'^e ~ 10 and u = 1/3 FQH system with 
a square unit cell. The horizontal axis is plotted in the Log 
scale to show the transition at low anisotropy. The system is 
in a Laughhn-liquid-like state for Ac < 2.0 fSec. IIII"^ . and 
becomes a quasi-lD crystal for Ac > 33.0 (Sec. IIII"C| . whose 
excitation gap scales as 1/^/34^. The intermediate regime is 
related to a Hall-smectic-hke phase, which is unstable towards 
the quasi-lD crystal and will be discussed in Sec. IIIIBl The 
inset shows a linear size-scaling for the minimum gap in the 
intermediate regime. 



in the inset reveals that the minimum gap can close in 
the thermodynamic limit, suggesting there might exist a 
switch between different order parameters ruling the sys- 
tem. Beyond Ac = 33.0 the energy gap decreases roughly 
as \j \fA'c^ indicating the regime of the quasi-lD repul- 
sion limit. We have studied other system sizes and found 
that these boundaries are size-dependent. But, in gen- 
eral, the ground state of the system maintains a three- 
fold degeneracy. This adiabatic transition with complex 
regimes typically occurs between distinct phases with a 
competition in the intermediate region. In the following 
subsections, we will focus on these different regimes in 
Ac and reveal an interesting competition between liquid 
and crystal phases. 



A. Anisotropic Laughlin Liquid at Small 
Interaction Anisotropy 

The energy gap plot suggests that the ground state at 
small interaction anisotropy is an incompressible liquid 
similar to the isotropic Laughlin state. The anisotropy in 
interaction, however, is expected to be imprinted, e.g., in 
the static structure factor of the resulting incompressible 
liquid. The projected static structure factor is defined 
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FIG. 2: (Color online) (a) 3D and contour plots of the struc- 
ture factor, (b) structure factor along q^ and qy axes, and (c) 
excitation spectrum along g^, and qy ajces for the N^ = 10 
system with Coulomb interaction anisotropy A^ — 1.8. To 
overcome the discrete momentum limitation, we use unit cells 
with different aspect ratios Ra = Lx/Ly in (b) to obtain more 
data points and a unit cell with Ra = 0.5 in (c). 



where |0) is the calculated ground state and r^ is the 
coordinate of the ith particle. 

In Fig. EJa), we draw the three-dimensional (3D) and 
contour plots of the structure factor for a calculated FQH 
state at Ac = 1.8. It exhibits a crater-like feature, which 
is similar to that of the isotropic liquid. However, the 
overall shape of the crater is deformed, stretching along 
the hard axis direction. Therefore, the elliptical symme- 
try replaces the circular symmetry in the isotropic liquid 
case. 

This anisotropic signature is more prominent in the 2D 
cuts along the two principal axes as shown in Fig. [2][b). 
We note that the structure factor behaves asymptotically 
as ^ q^ in the long wavelength limit. This agrees with 
the single mode approximationiS. (SMA) for incompress- 
ible liquid. However, the prefactor of the quartic term is 
orientation dependent, revealing the anisotropic nature 
of the structure factor. According to Ref. [y], the ratio of 
prefactors at Qx and Qy axes is equal to (l/^^)"*, where 
the parameter A\ defines an intrinsic metric, describ- 
ing how the correlated quasi-particles bind to each other 
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FIG. 3: (Color online) (a) Pair correlation function along x 
and y axes for the N^ ~ 10 system at interaction anisotropy 
Ac — 1.8. The dotted-dashed line represents the correlation 
function at j4c = 1 for comparison. The inset reveals the ar^ 
behavior at small r for the anisotropic case, (b) A linear fit 
for the prefactor ^/a at small anisotropy. 



in the anisotropic environment. The fitting lines in the 
plot have revealed A*j^ ^ 1.45 for interaction anisotropy 
Ac = 1.8. The peaks in the orientation-dependent plots 
represent the crater ridge in Fig. [2lja). 

According to the SMA, the maximum in the struc- 
ture factor corresponds to a minimum gap in the excita- 
tion spectrum, or the roton minimum, which corresponds 
to the excitonic binding of the neutral quasiparticle- 
quasihole pairs. ^^ This is evident in Fig. [He), where 
we plot the orientation-dependent low-energy excitation 
spectra in the momentum space. The location of the ro- 
ton minimum is sensitive to the direction, but the gap 
value is less sensitive. The ratio of the roton-minimum 
locations along q^ and Qy axes is found close to A*j^ as ex- 
pected and these two locations match the peak locations 
of the structure factor in Fig. [DJb). 

The isotropic Laughlin wavefunction, with order-3 ze- 
ros at the locations of other particles, triumphed in 
the explanation of the isotropic incompressible liquids 
of 1/ = 1/3 FQH system. Corresponding to the deformed 
electron-hole correlation from the anisotropic interaction, 
order-3 zeros in the wavefunction are expected to split. 
Several works have suggested that the relative coordinate 
part of the anisotropic wavefunction^iiSr— has the form 



w(zj) -- ^(^j - Zj){{^l 



? + 4) 



(5) 



i<j 



where Zi is the complex coordinate of the ith particle and 
zq is a complex constant related to the splitting of the 
zeros due to anisotropy. This zero-splitting effect in the 
wavefunction can be detected using the pair correlation 
function defined asii 



5(i-) 



Lx-Ly 



Nc{Nc - 1) 



\Y,6ir-ir,-v,))\0). (6) 



In Fig. Elja), we plot the pair correlation functions 
along X and y directions for the ly = 1/3 FQH state with 
the interaction anisotropy Ac — 1.8. The two curves 
are distinguishable from their isotropic counterparts. We 
point out that the curves behave asymptotically as ar^ 
in the limit of r — )> 0, with the prefactor a oc |zo|^- This 
is entirely different from the isotropic Laughlin wave- 
function, which exhibits a r^ asymptotic behavior in its 
pair correlation function. The nonmonotonic behavior in 
g{y) at small r region, which manifests itself more clearly 
at a larger Ac, is also consistent with the zero-splitting 
scenarios. i^ii^ We point out that for a suitable deformed 
model wavefunction, there is also an additional contribu- 
tion to the Gaussian Landau level form factor, ^'^* which 
can be observable in disk geometry with a boundary. 

In Fig. El^b) , we plot the square root of the prefactor a 
at several small anisotropy. The linear fit of y/a to Ac is 
expected as {Ac — 1) [or {y/A^~ 1)], which characterizes 
the perturbation away from the isotropic point. However, 
the resulting nonzero intercept at Ac = I suggests that 
we may have overestimated the prefactor, possibly due 
to the higher order contributions at small Ac- 
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FIG. 4: (Color online) (a)-(c): Structure factors for the 
ground states of the A^e = 10 system at Ac = 5.0, 11.0, and 
20.0, respectively. Note that the locations of the twin peaks 
for Ac = 5.0 (CDW-like peaks) differ from those for Ac = 20.0 
(Crystal peaks). Both types of peaks coexist at Ac — 11.0. 
(d) Peak values of the CDW-like and Crystal peaks are plot- 
ted as a function of the interaction anisotropy in a range of 
l<Ac< 50. 



B. Hall-Smectic-like Phase in the Intermediate 
Interaction Anisotropy Regime 

The explicit construction^, of the model wavefunction 
by unimodular transformation on disk geometry suggests 
that the geometrical description of the quantum Hall sys- 
tem accepts the following deformation of the isotropic 
Laughlin state (i.e., 7 = 0) 



*-n 



i<j 



where z,- 



127* 



I-I7P 



-i:.7^?/4p-E. 



V4 



(7) 



Zi — Zj and 7 characterizes the amount of 
mixing between the guiding center creation and annihi- 
lation operators in the unimodular transformation. Note 
that the model wavefunction is expected to be valid for 
small 7. In the present parametrization 7 — \/A^ — 1 is 
real. We can postulate the breakdown criterion for the 
anisotropic Laughlin liquid to be 



127* 



2tt 



(8) 



i.e., the area occupied by a set of three splitting zeros 
is the average area per particle. This suggests that at 
the breakdown A*^ k, 2 (i.e., Aj, « 3 according to the 
estimation in the subsection IIII D[) , consistent with the 
onset of the rapid decrease of the excitation gap. 

In other words, the anisotropic Laughlin liquid is stable 
when the long-distance (i.e., at average particle spacing) 



behavior of the Jastrow factor is still as z. 
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excitations, which becomes anisotropic. When the liquid 
phase breaks down, it cannot sustain further anisotropy 
by the spatial deformation in the roton spectrum. 

One possible outcome of the system after this break- 
down is that the mode at the roton minimum goes softer, 
developing some charge-density-wave (CDW) order. Due 
to the orientation effect of the anisotropy, this CDW 
is expected to be unidirectional (stripe-like) and the 
characterizing sharp peaks in the structure factor are 
along the stretching direction. This is clearly visible 
in the structure factor at Ac — 5.0 in Fig. lU^a). The 
background in the structure factor resembles that of an 
anisotropic Laughlin liquid, but its peak value is signif- 
icantly smaller than the two sharp peaks along the q^ 
axis. The CDW twin peaks correspond to a period in real 
space, which can be roughly anticipated as the splitting 
of zeros |zo| = \/6 at the critical A*j^ = 2. The peak value 
(subtracting background) in the structure factor sufhces 
as the order parameter. The plot in Fig. IH^d) shows that 
this CDW-like order parameter rules the system in the 
regime 2.0 < Ac < 8.0. 

We term this phase, which breaks one-dimensional 
translational symmetry, as a Hall-smectic-like phase since 
we speculate that it is related to the Hall smectic dis- 
cussed earlier in the context of liquid crystal phases 
in the FQH system.— "—li^"— The rise of the smectic 
phase softens the magnetoroton mode and appears to 
be responsible for the reduction of the excitation gap for 
2.0 < Ac < 8.0 as shown in Fig. [TJ As discussed in 



five excitation of the liquid is the neutral magnetoroton 



Ref. [ll| the transition from the Laughlin liquid to the 
Hall smectic can be second order and its critical behavior 
is in the XY universality class. 

Beyond Ac — 8.0, the reverse trend in the excitation 



gap as a function of Ac indicates that the system is un- 
der the influence of a distinct mechanism. This is evi- 
dently shown in the structure factor plot of Fig. IDJb) at 
Ac — 11.0. Two additional peaks along the q^ axis are 
clearly visible with the different wavevectors from the 
CDW-like twin peaks. These additional twin peaks are 
corresponding to the unidirectional crystal order in the 
quasi- ID repulsion limit that we will discuss in subsection 
unci The peak value of them is plotted as the crystal or- 
der parameter in Fig. IDJd). There we can see that the 
crystal order parameter is continuously increasing with 
the interaction anisotropy. The crossover for the compe- 
tition with the CDW-like order occurs around Ac = 13.0. 
For larger anisotropy the crystal order dominates as illus- 
trated in Fig. mjc) at Ac = 20.0, where only the crystal 
peaks remain. Thus, the Hall-smectic-like phase is found 
unstable towards a ID crystal. 



C. Quasi-ID Crystal in the Large Anisotropy Limit 

The crystal phase at the large anisotropy limit can be 
probed using the charge distribution. In Fig. EJa), we 
plot the average LL orbital occupations and the charge 
density along the hard axis at ^c = 40.0. The charge den- 
sity appears smoother as an integral from local Gaussian 
wave packets over orbital occupations and a 2% fluctu- 
ation above the background can be observed in the ex- 
aggerated plot. Both the charge occupation and density 
fluctuate along the hard axis with the crystalline period 
A* = Lx/Nc- The maximum charge occupation is close 
to unity as expected in the ultimate ID crystal limit. The 
2D distribution of the charge density reveals that the sys- 
tem is a unidirectional crystal with each electron spread- 
ing into a stripe perpendicular to the hard axis. For basis 
states in the torus geometry, the guiding-center coordi- 
nate along the hard axis is coupled with the momentum 
along qy axis." Thus, the calculated ground states are 
expected to carry a period of 2t: /Xy in the momentum 
space along the qy axis, where Aj, = Ly/i. The structure 
factor plot in Fig. Eljb) demonstrates this characteristic 
order along qy axis. As its real space counterpart, the 
pair correlation plot in Fig. [SJc) shows oscillation in y 
direction with a period Aj,. 

The above results support that at filling v = 1/3 the 
ground state of the system is a crystal in the large inter- 
action anisotropy limit and the system undergoes some 
transition from an incompressible liquid to a solid as 
anisotropy increases. A similar story has been discussed 
in an isotropic FQH system with extreme geometry, such 
as in a thin torus or a cylinder limit f22i"— and in a recent 
work25, on the graphene ribbon with flat bands. They 
can be explained under the same principle in Ref. [2J| 
by sorting Hamiltonian. When the interaction anisotropy 
increases, the repulsion-related diagonal terms dominate, 
which has the similar effect as geometry on the isotropic 
FQH and as the local orbital expansion on the flat-band 
graphene ribbon. The low-energy physics is governed by 
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FIG. 5: (Color online) For ground states of A^e = 10 system 
at Ac = 40.0: (a) Average occupations (scattered squares) 
of the Landau-level orbital at guiding centers and the charge 
density (solid line) are plotted along the hard axis. A period 
of A* = Lx/Ne is clearly visible. The density values have been 
shifted and exaggerated to emphasize a ~ 2% fluctuation over 
the background, (b) Structure factor with sharp twin peaks 
in (fa; direction and periodic osciUation in qy direction, (c) 
Pair correlation function is plotted along two principal axes. 
The function g{y) shows a osciUation with the period Ly/3. 



the strong repulsion so that the system tends to form 
crystal. At small anisotropy, the hopping-related off- 
diagonal terms are comparable and screen the repulsion, 
resulting in the liquid phase. 



D. Generalized Variational Laughlin State 

In the discussion above, we have seen that the isotropic 
Laughlin wavefunction is insufficient to fully capture the 
features of an anisotropic FQH system. For such a system 
at the lowest LL filling i/ = l/q, Haldane has suggested 
to use a family of Laughlin states,^ which is generally de- 
fined as the densest zero-energy eigenstate of a projected 
two-body anisotropic Hamiltonian: 



H,{Al) - Y. P^i^L). 



(9) 



m<q 



For the fermion system with an odd denominator g, m 
are limited to be odd. This Hamiltonian is a truncated 
summation over anisotropic pair interactions: 






(10) 



for two particles with the relative angular momentum 
of mh in the guiding-center coordinates. In the above 
expression, Lra{x) are mth Laguerre polynomials and 
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FIG. 6: (Color online) Wavefunction square overlap (lines) 
and expected value of Coulomb energy (symbols) as a function 
of the variational parameter Al for v — 1/3 FQH system 
at interaction anisotropy Ac — 1.8. The energy values have 
been shifted and enlarged to emphasize that the location of 
the minimum coincides with that of the largest overlap. The 
comparison of different sizes of N,; = 6, 8, 10, and 12 shows a 
weak size dependence only. 



Q(Al) = Jql/AL + AlqI, which, like in Eq. ([3]), de- 
fines a wavefunction metric parameterized by Aj^. The 
parameterized Laughlin states $(Al) satisfy 



P,„(Ai)|^(^i))=0,TO<g. 



(11) 



The isotropic Laughlin wavefunction corresponds to the 
Laughlin state with Al = 1. With this family of pa- 
rameterized states, we are able to variationally approxi- 
mate the ground state of a FQH system with anisotropic 
interaction. According to Haldane's proposal, if the 
mass or orbital metric (in our case, an identity ma- 
trix for isotropic mass) is different from the interaction 
metric (parameterized by Ac), the resulting variational 
state ^[A*i) should be described by a metric interpo- 
lating the mass metric and the interaction metric, i.e., 
1 < A*]^ < Ac- This intrinsic metric describes how corre- 
lated quasi-particles effectively feel each other in a such 
anisotropic FQH system. 

In Fig. [HI we study the anisotropic i^ = 1/3 FQH sys- 
tem with the Coulomb anisotropy Ac — 1.8. The optimal 
Laughlin state ^(yl^) is obtained by tracing either the 
maximum of the wavefunction overlap or the minimum 
of the expected Coulomb energy 

Ea{AL) = (vI/(Ai)|i7c(Ac)|*(Ai)). (12) 

The optimal parameter is found at A'^ ~ 1.43, which 
is weakly size-dependent. This parameter is indeed 
an intermediate value between unity and the Coulomb 
anisotropy as expected.^ It also agrees with the intrinsic 
metrics through the analysis of the anisotropic structure 
factor in the subsection IIII Al The overlaps between the 
optimal Laughlin state and the exact ground states are 



FIG. 7: (Color online) For a Coulomb system with interac- 
tion anisotropy Ac = 1.8: (a) Effective pseudopotentials Vm 
with a given parameter Al = 1.47 is plotted against m. (b) 
The pseudopotential Vi is plotted as a function of Al- The 
maximum value occurs around Al = 1.47. 



larger than 90% for various system sizes, which supports 
the validity of the variational state. We also note that 
the expected Coulomb energy quadratically approaches 
its minimum, which suggests a linear approximation of 
the anisotropic Laughlin state with Al in the liquid phase 
regime. 

To gain a further understanding for the validity of 
this variational approach, we approximately expand the 
Coulomb interaction in the anisotropic pair interactions 
as: 



Hc{Ac) « J2 ^™(^-' AL)PmiAL), 



(13) 



where the average expansion coefficients Vm define the 
effective anisotropic pseudopotentials in a form of 



14,,= 



d9 



dx 



Lm{x^)e 



-F{e,x)/2 



27T^/G{0) 



with 



and 



Gi9) = (AJAl) cos^ 9 + {Al/Ac) sin^ 9 



F{9,x) =x^{l + Al cos^ 9 + sin^ 9/ Al). 



(14) 



(15) 



(16) 



An example of Vm is plotted in Fig. [Tja) with Ac — 
1.8 and Al = 1.47. These pseudopotential parameters 
are found positive and monotonously decrease with m, 
consistent with the long range behavior of the Coulomb 
repulsion. 

Given the Coulomb interaction anisotropy Ac , in prin- 
ciple we can have a family of pseudopotential sets pa- 
rameterized by {^l}, which are associated with different 
values. The set of Vm with the maximum pseudopoten- 
tial values up to the mth order is most promising for the 
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FIG. 8: (Color online) Contour plot of wavefunction square 
overlap between the Coulomb ground state (parameterized by 
Ac) and the generalized Laughlin state (parameterized by Al) 
for the A^e = 10 system. Each dot represents the maximum 
of the wavefunction overlap for a given Ac- The solid curve 
stands for the estimated fit A^ as a function of Ac using the 
maximum Vi criterion. 



truncation-based variational approach to work. Thus, for 
the 1/ = 1/3 system, we could use the maximum of Vi as 
a criterion to examine the optimal parameter of the vari- 
ational Laughlin state, A^. This size-independent condi- 
tion 



ViiAl) > Vi{Ai 



(17) 



serves as a semi-analytic estimation to the intrinsic geom- 
etry parameter A^ of the z^ = 1/3 anisotropy system. As 
shown in Fig. EJ^b), the estimated optimal parameter is 
around A"^ — 1.47, matching the value found previously 
through the finite-size calculation. 

In the 2D contour plot of Fig. [SJ we show a compre- 
hensive map of the wavefunction overlap as a function 
of both parameters Ac and Al ■ The isotropic Laughlin 
wavefunction has a continuously decreasing overlap with 
the calculated FQH state when the Coulomb anisotropy 
increases, indicating again its insufficiency in describing 
the anisotropic FQH system. Instead, We note that at 
the optimal parameters, the local maximum of the wave- 
function overlap is larger than 90% in the full range of 
the Coulomb anisotropy parameter {even when the sys- 
tem is not in the liquid phase). This justifies that the 
family of variational Laughlin states, which have larger 
wavefunction overlaps, could serve a better description of 
the anisotropically interacting FQH system. We also no- 
tice that in the vicinity oi Ac — 1, the optimal parameter 
of the Laughlin state is linearly related to the interaction 
anisotropy. As a comparison, we plot the estimated op- 
timal parameters A^ from Eq. PT)) as a function of Ac- 
We find they agree well with the finite-size results in the 
liquid and crystal phases. But in the intermediate regime 



the results show some deviation. This suggests a more 
careful handling beyond the simple criterion of Eq. p7)) 
is needed for the intermediate region. 



IV. SUMMARY 

In conclusion, we have studied the effect of anisotropic 
Coulomb interaction on the ground state of the 1/3 fill- 
ing fractional quantum Hall system. We find that at 
weak anisotropy the Laughlin state remains to be a valid 
description, although the structure factor and pair corre- 
lation function exhibit anisotropy. Our calculations sup- 
port the recent proposal of Haldane on the geometric 
description of the fractional quantum Hall state. ^ In par- 
ticular, the order-3 zeros in the wavefunction split into 
three distinct zeros with a splitting distance related to the 
anisotropyi^ We have compared the ground state wave- 
function of the anisotropic Coulomb interaction with a 
family of single-parameter variational Laughlin states. 
The latter are obtained by deforming the projection (i.e., 
Vi only) Hamiltonian for the isotropic Laughlin state. 
We have determined the variational parameter by min- 
imizing the variational ground state energy or by max- 
imizing the wavefunction overlap. In addition, we also 
propose an effective analysis to estimate the optimal 
variational parameter. The liquid phase breaks down 
when anisotropy increases and a Hall-smectic-like order 
emerges. Finally, at strong interaction anisotropy, which 
is more of theoretical interest, the ground state exhibits 
a compressible one-dimensional crystal phase. Interest- 
ingly, the ground state obtained by solving deformed pro- 
jection Hamiltonian remains to be a good description 
(overlap greater than 90%) throughout the liquid to solid 
transition. 

The anisotropy in Coulomb interaction provides a new 
route in probing the intrinsic metric of the fractional 
quantum Hall state in its geometrical description, as 
pointed out by Haldane.i In the present study we have 
revealed that the intrinsic metric or the wavefunction 
anisotropy are indeed different from the Coulomb metric 
or the dielectric tensor anisotropy, although they appear 
to be linearly proportional to each other in the vicinity of 
the isotropic Laughlin state. The linearity we have found 
supports the proposal to use a single-parameter to con- 
struct unimodularly deformed wavefunctions to describe 
the effect of interaction anisotropy in disk geometry, - In 
fact, the family of deformed wavefunctions contain the 
same anisotropic Jastrow factor in their relative coor- 
dinate part,^*^ which explicitly splits the order-3 zeros 
in the Laughlin liquid. Therefore, the emergence of the 
Hall smectic phase, as analyzed in the effective field the- 
orieSfii^ at larger anisotropy is also a support of the 
geometrical description of FQH states. 

Near the completion of this work, we note a very recent 
work2£ which discusses the anisotropic FQH system with 
the anisotropic band mass. There, the persistent energy 
gap and anisotropic roton-minimum excitation at small 
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anisotropy parameters are reported, consistent with our 
results in the hquid phase. 
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